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Abstract. Current status and future prospects of the structure functions and parton distribution 
studies are presented. 

INTRODUCTION 

The study of nucleoli's structure functions and parton distributions is an active area 
of research in nuclear and particle physics. The parton distributions address both the 
perturbative and nonperturbative aspects of QCD, and they also provide an essential 
input for describing hard processes in high-energy hadron collisions. As a result of 
several decades 's intense effort, the unpolarized proton structure functions have been 
well mapped out over a broad range of Q 2 and Bjorken-x While these data are invaluable 
for testing QCD and for extracting various parton distributions, several questions remain 
unanswered. For example, the unexpected finding of the flavor asymmetry of the light- 
quark sea (u, d) suggests that other aspects of the flavor structure, such as possible 
asymmetry between the s and s sea quark distributions and the bahavior of valence d/u 
ratio at large x, need to be examined. The issue of quark-hadron duality, as reflected 
in the intriguing similarity between the structure functions measured at the resonance 
region and at the DIS region, also requires further studies. 

Remarkable progress in the study of spin-dependent structure functions has been 
made since the discovery of the "proton spin puzzle" in the late 1980's. Very active 
spin-physics programs have been pursued at many facilities including SLAC, CERN, 
HERA, JLab, and RHIC. The polarized DIS data now cover a sufficiently broad Q 2 
range for scaling-violation to be observed. In recent years, new experimental tools such 
as semi-inclusive polarized DIS, polarized proton-proton collision, and deeply exclusive 
reactions have been employed to address the major unresolved question in spin physics: 
How is the proton's spin distributed among its various constituents? 

On the theory front, the formulation of the generalized parton distributions as well as 
the identification of various kj (intrinsic transverse momentum of partons)-dependent 
structure and fragmentation functions have opened exciting new directions of research. 
Furthermore, important progress in the Lattice calculations for the moments of various 
parton distributions and in the extrapolations to their chiral limits has been made. 

In this review I will focus on recent progress in the following areas: 

• Flavor structure of parton distributions 

• Transition from high-Q 2 to low-<2 2 

• Novel distribution and fragmentation functions; Generalized parton distributions 
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FIGURE 1. Left panel: Cross section ratios of p + d over 2(p + p) for Drell-Yan, J/^f*, and T production 
from FNAL E866. Right panel: Comparison of E866 d — u data with calculations from various models 0]. 

FLAVOR STRUCTURES OF PARTON DISTRIBUTIONS 

d/u flavor asymmetry 

The earliest parton models assumed that the proton sea was flavor symmetric, even 
though the valence quark distributions are clearly flavor asymmetric. The flavor sym- 
metry assumption was not based on any known physics, and it remained to be tested. 
Under the assumption of a u, d flavor- symmetric sea in the nucleon, the Gottfried Sum 
Rule [1], Iq = /J [F%(x, Q 2 ) -F£(x, Q 2 )) /xdx=\ /3, is obtained. The NMC collabora- 
tion determined the Gottfried integral to be 0.235 ± 0.026, significantly below 1/3. This 
surprising result can be explained by a large flavor asymmetry between the u and the d. 

The x dependence of dj u asymmetry has been determined by proton-induced Drell- 
Yan (DY) as well as semi-inclusive DIS measurements. Figure 1 shows that the Fermilab 
E866 [2] DY cross section per nucleon for p + d clearly exceeds p + p, and it indicates 
an excess of d with respect to u over an appreciable range in x. In contrast, the o(p + 
d)/2o(p + p) ratios for J and T production, also shown in Fig. 1, are very close 
to unity. This reflects the dominance of gluon-gluon fusion process for quarkonium 
production and the expectation that the gluon distributions in the proton and in the 
neutron are identical. 

Many theoretical models, including meson cloud model, chiral-quark model, Pauli- 
blocking model, instanton model, chiral-quark soliton model, and statistical model, have 
been proposed to explain the d/u asymmetry. For recent reviews, see O, |4|]. These 
models can describe the d — u data very well, as shown in Fig. 1. However, they all 
have difficulties explaining the d/ u data at large x (x > 0.2). The new 120 GeV Fermilab 
Main Injector and the proposed 50 GeV Japanese Hadron Facility present opportunities 
for extending the dj u measurement to larger x (0.25 < x < 0.7). 



Models in which virtual mesons are admitted as degrees of freedom have implications 
that extend beyond the d,u flavor asymmetry addressed above. They create hidden 
strangeness in the nucleon via such virtual processes as p — > A + K + ,E + K, etc. Such 
processes are of considerable interest as they imply different s and sparton distributions 
in the nucleon, a feature not found in gluonic production of ss pairs. 

A difference between the s and s distribution can be made manifest by direct mea- 
surements of the s and s parton distribution functions in neutrino DIS. A fit to the CDHS 
neutrino charged-current inclusive data together with charged lepton DIS data found ev- 
idence for Jq s(x)dx > Jq s{x)dx o]. However, an analysis 10] of the recent CCFR and 
NuTeV v(v)N — > n + /l~x dimuon production data 0] favored J l s(x)dx < J l s(x)dx 

(Jq (s(x) — s(x))dx = —0.0027 ±0.0013). To better determine the s/s asymmetry, an 
NLO analysis is currently underway [8]. Violation of the s/s symmetry would have im- 
pact on the recent extraction [9] of s'm 2 6w from the CCFR/NuTeV vN scattering data. 

Asymmetry in the s,s distributions can also be revealed in the measurements of the 
strange quark's contribution to the nucleon's electromagnetic and axial form factors. 
These "strange" form factors can be measured in neutrino elastic scattering [10] from 
the nucleon, or by selecting the parity-violating component of electron-nucleon elastic 
scattering. Two completed parity- violating experiments 111 ill 1211 suggest small contribu- 
tions of strange quarks to nucleon form factors. Several new experiments are underway 
at JLab and MAMI to measure parity- violating asymmetry at various kinematic regions. 



Flavor structure of polarized nucleon sea 

The flavor structure and the spin structure of the nucleon sea are closely connected. 
Many theoretical models originally proposed to explain the dj u flavor asymmetry also 
have specific implications for the spin structure of the nucleon sea. In the meson- 
cloud model, for example, a quark would undergo a spin flip upon an emission of a 
pseudoscalar meson (w — > %°{uu,dd) + u^, — > n + (ud) + d\ u) — > K + +s+, etc.). 
The antiquarks (u, d, s) are unpolarized (Au = Ad = As = 0) since they reside in spin-0 
mesons. The strange quarks (s), on the other hand, would have a negative polarization. 

In the chiral-quark soliton model lfl3l fbUl . the polarized isovector distributions 
Aii(x) — Ad(x) appears in leading-order (N^) in a l/N c expansion, while the un- 
polarized isovector distributions u(x) — d(x) appear in next-to-leading order (N c ). 
Therefore, this model predicts a large flavor asymmetry for the polarized sea 
[Au(x) —Ad(x)) > [d(x] —u(x)\. 

The HERMES collaboration has recently reported the extraction of Aii(x), Ad(x), 
and As(x)(— As(x)) using polarized semi-inclusive DIS (SIDIS) data lfl5ll . Although 
the statistics are still limited, the HERMES results for Au, Ad, Au — Ad, as shown in 
Fig. 2, are all consistent with being zero. In particular, there is no evidence for a 
large positive Au(x) — Ad(x) asymmetry as was predicted [16] by the chiral quark 
soliton model. Figure 2 also shows that As tends to be positive, in contrast to the 
predictions of a negative polarization of the strange sea in the analysis of inclusive DIS 
and hyperon decay data assuming SU(3) symmetry. However, the HERMES result of 
As = 0.03 ± 0.03 ± 0.01 over 0.023 < x < 0.3 is not in disagreement with the inclusive 




FIGURE 2. Quark and antiquark polarizations extracted from the HERMES SIDIS data Pi . 

DIS result of (As + As)/2 ~ -0.02 0. 

Another promising technique for measuring sea-quark polarization is W -boson pro- 
duction 111 811 at RHIC. The longitudinal single- spin asymmetry for W production in po- 
larized p + p — > W ± + x gives a direct measure of sea-quark polarization. The RHIC 
^-production and the HERMES SIDIS measurements are clearly complementary tools 
for determining polarized sea quark distributions. 



d/u ratio at large x 

Another quantity related to the flavor symmetry of the proton is the d/u ratio at large 
x. Assuming SU(2) sp j>, x SU(2)// avor symmetry, the proton wave function is given as 

\p >T = -t=" T Ms=o,s z =o + — ]=u t (ud)s=i,s z =o ~ j« i {ud) s =\,s z =\ 

V2 vlo J 

1 y/2 

--d t {uu)s=i,Sz=0 + —d I {uu)s=\,s z =\ (1) 

The neutron wave function is readily obtained from u <-> d interchange. In nature, the 
SU(2). sp ,„ x SU(2) fi avor symmetry is clearly broken, as evidenced by the large N — A 
mass splitting. The dynamic origins of this symmetry breaking remains unclear. Close 
and Carlitz I19L 12011 argued that the dominance of the 5 = diquark configuration over 
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FIGURE 3. Left panel: u and d valence quark densities obtained from HI charged-current measure- 
ments 1 25]. Right panel: Comparison of the E866 pp and pd DY cross sections with PDF calculations [26]. 



the 5=1 configuration would account for the N — A mass splitting as well as the 
SU(2) x SU(2) symmetry breaking. An alternative suggestion, based on perturbative 
QCD, was offered by Farrar and Jackson 12 ill . They pointed out that the spin-aligned 
diquark configuration with Sz = 1 is suppressed since only longitudinal gluons can be 
exchanged. A similar result was also obtained by Brodsky et al. 12211 using counting 
rule argument. It is straightforward to show that in the x — ► 1 limit, the different models 
predict the folllowing values for various ratios: 
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The distinct predictions for F^/FJ? from various models could be tested against DIS 
experiments. However, there exist considerable uncertainties in the extraction of F£ 
from the measurement of F^ ■ Depending on the treatment of the nuclear effects in the 
deuteron, very different values for F% /Ff (and d/u) were obtained at large x 112311 . It 
is clearly desirable to measure d/u without the need to model nuclear effects in the 
deuteron. One method is to measure the charge asymmetry of W production in p — p 
collision. Indeed, the CDF data 12411 on the W charge asymmetry have already provided 
useful constraints on the d/u ratio. 

The d/u ratio can also be probed by measuring the e~ p — > v e x and e + p — > V e x 
charged-current DIS, where the underlying processes are e~u — > v e d and e + d — > v e u, 
respectively. The recent HI charged-current data 12511 . shown in Fig. 3, indicate that 
the u quark density at large x(x = 0.65) is smaller than expected from the current 



PDF parametrization. Very recently, the Fermilab E866/NuSea collaboration reported 
the absolute Drell-Yan cross sections of 800 GeV p + p and p + d 112611 . As shown in 
Fig. 3, the data fall below the PDF predictions at large x (up to x = 0.8). The HI and the 
E866 results suggest that u quark density at large x might be smaller than expected from 
current PDFs. This clearly would impact on the d/u ratio at large x, as shown in a recent 
global PDF analysis [& 

The uncertainties involved in the extraction of F 2 " from Fj data can be greatly reduced 
using the technique of neutron-tagging. A new experiment [28] has been proposed at the 
JLab Hall-B to detect e~d^ e~ px, where a low-energy recoiled proton will be measured 
in coincidence with the (e, e') scattering. Using this method, the F^ /F 2 P rat ^° over me 
range 0.2 < x < 0.7 could be determined with small systematic uncertainties. 

TRANSITION FROM HIGH-0 2 TO LOW-g 2 
Quark-hadron duality 

The recent studies at JLab of the spin-averaged and spin-dependent structure functions 
at low Q 2 region have shed new light on the subject of quark-hadron duality. Thirty years 
ago, Bloom and Gilman [29] noticed that the structure functions obtained from deep- 
inelastic scattering experiments, where the substructures of the nucleon are probed, are 
very similar to the averaged structure functions measured at lower energy, where effects 
of nucleon resonances dominate. This surprising similarity between the resonance elec- 
troproduction and the deep inelastic scattering suggests a common origin for these two 
phenomena, called local duality. 

Recently, high precision data ll3(ill from JLab have verified the quark-hadron duality 
for spin-averaged scattering on proton and deuteron targets. For Q 2 as low as 0.5 GeV 2 , 
the resonance data are within 10% of the DIS results. When the mean F2 curve from 
the resonance data is plotted as a function of the Nachtmann variable, £ = 2x/(l + 
■s/l + 4M 2 x 2 / Q 2 ) , it resembles the xF$ structure function obtained in neutrino scattering 
experiments. Since xF^ is a measure of the valence quark distributions, this suggests that 
the F2 structure function at low Q 2 originates from valence quarks only. 

The study of quark-hadron duality was recently extended to other structure functions. 
Results from HERMES 113 ill show that duality is also observed for the spin-dependent 
quantity Aj\ Another recent result from JLab shows that the nuclear modifications to 
the unpolarized structure functions in the resonance region are in surprisingly good 
agreement with those measured in DIS 13211 . 

I~i (Q 2 ) at low Q 2 and the generalized GDH integral 

The extensive data on gi(x, Q 2 ) allow accurate determinations of the integrals 
rf n (Q 2 ) = Jq g^ n {x, Q 2 )dx for the proton and the neutron, as well as T\(Q 2 ) -T'l(Q 2 ). 
While the values of and T" are different from the predictions of Ellis and Jaffe who 




assumed SU(3) flavor symmetry and an unpolarized strange sea, the data are in good 
agreement with the prediction of the Bjorken sum rule. 

How does Ti(Q 2 ) evolve as Q 2 — > 0? This question is closely related to the 
Gerasimov-Drell-Hearn (GDH) sum rule: 

[<Tl/2(v)-a 3/2 (v)] — = --^ r K 1 . (2) 

The GDH sum rule, based on general physics principles (causality, unitarity, Lorentz and 
gauge invariances) and dispersion relation, relates the total absorption cross sections of 
circularly polarized photons on longitudinally polarized nucleons to the static properties 
of the nucleons. In Eq.El O1/2 and a 3 / 2 are the photo-nucleon absorption cross sections 
of total helicity of 1/2 and 3/2, v is the photon energy and Vo is the pion production 
threshold, M is the nucleon mass and K is the nucleon anomalous magnetic moment. 
The GDH integral in Eq. |2] can be generalized from real photon absorption to virtual 
photon absorption with non-zero Q 2 : 



igdh(q 2 ) - /V 1/2 (v,e 2 ) -<r 3/2 (v,e 2 )]^ = ^^(q 2 ) 



(3) 



Eq. [3] shows that the <2 2 -dependence of the generalized GDH integral is directly related 
to the <2 2 -dependence of Ti. The GDH sum rule (Eq.EJ) predicts = at Q 2 = with 
a negative slope for dr^(Q 2 )/dQ 2 and is known to be positive at high Q 2 , therefore, 
rf(<2 2 ) must become negative at low Q 2 . 

The GDH integrals at low Q 2 have recently been measured in several experiments 



at JLab ||3J,|33|] and HERMES [35]. Results from a JLab Hall-B measurement ||33j] 
of T^(Q 2 ) are shown in Fig. 4. These data indeed show that T 1 ^ changes sign around 
Q 2 = 0.3 GeV 2 . The origin of the sign-change can be attributed to the competition 



between A(1232) and higher nucleon resonances. At the lowest Q 2 , the A(1232) has 
a dominant negative contribution to Ff. However, at larger Q 2 , higher mass nucleon 
resonances take over to have a net positive T^. 

Results [34] from a JLab Hall-A measurement of the generalized GDH integral for 
neutron using a polarized 3 He target are shown in Fig. 4. In contrast to the proton case, 
the strong negative contribution to the GDH integral from the A(1232) resonance now 
dominates the entire measured Q 2 range. Future experiments at JLab will extend the 
measurements down to Q 2 = 0.02 GeV 2 in order to map out the low Q 2 behavior of the 
neutron and proton generalized GDH integrals. 



NOVEL DISTRIBUTION AND FRAGMENTATION FUNCTIONS 

In addition to the unpolarized and polarized quark distributions, q(x, Q 2 ) and Aq(x, Q 2 ), 
a third quark distribution, called transversity, is the remaining twist-2 distribution yet to 
be measured. This helicity-flip quark distribution, dq(x, Q 2 ), can be described in quark- 
parton model as the net transverse polarization of quarks in a transversely polarized 
nucleon. Due to the chiral-odd nature of the transversity distribution, it can not be 
measured in inclusive DIS experiments. In order to measure 8q(x, Q 2 ), an additional 
chiral-odd object is required. For example, the double spin asymmetry, Att, for Drell- 
Yan cross section in transversely polarized pp collision, is sensitive to transversity since 
Att ~ Hi e f^^i( x i)^^i( x 2)- Such a measurement could be carried out at RHIC ifl^l . 
although the anticipated effect is small, on the order of 1 — 2%. 

Several other methods for measuring transversity have been proposed for semi- 
inclusive DIS. In particular, Collins suggested [36] that a chiral-odd fragmentation func- 
tion in conjunction with the chiral-odd transversity distribution would lead to a single- 
spin azimuthal asymmetry in semi-inclusive pion production. 

The HERMES collaboration recently reported IB7I1 observation of single-spin az- 
imuthal asymmetry for charged and neutral hadron electroproduction. Using unpolarized 
positron beam on a longitudinally polarized hydrogen and deuterium targets, the cross 
section was found to have a sin0 dependence correlating with the target spin direction. 
is the azimuthal angle between the pion and the (e, e') scattering plane. This Single- 
Spin-Asymmetries (SSA) can be expressed as the analyzing power in the sin0 moment, 
and the result is shown in Fig. 5. The sin0 moment for an unpolarized (U) positron 
scattered off a longitudinally (L) polarized target contains two main contributions 

1 ~y q q 

where Si and St are the longitudinal and transverse components of the target spin 
orientation with respect to the virtual photon direction. For the HERMES experiment 
with a longitudinally polarized target, the transverse component is nonzero with a mean 
value of St ~ 0.15. The observed azimuthal asymmetry could be a combined effect 
of the h\ transversity and the twist-3 hi distribution. Recently, another mechanism 
involving a chiral-even T-odd Sivers distribution function [38] was shown to contribute 
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FIGURE 5. Analyzing power in the sin0 moment from HERMES [ 3 



to azimuthal asymmetry |39[ 4(j]. For a longitudinally polarized target the Collins and 
the Sivers mechanisms can not be distinguished. 

If the azimuthal asymmetry observed by HERMES is indeed caused by the h\ 
transversity, a much larger asymmetry is expected for a transversely polarized target. The 
HERMES and COMPASS collaborations have collected polarized SIDIS using trans- 
versely polarized hydrogen and 6 LiD targets, respectively. These data would shed much 
light on the origins of the SSA and could also disentangle the Sivers effect from the 
Collins effect. The Collins effect has a sin(<^ + § l s ) dependence while the Sivers effect 
is proportional to sin(<^' — where 0j — (j) s — (j) 1 is the angle between target spin and 
the lepton scattering plane. For longitudinally polarized target 0j = and the two effects 
have identical dependence. For transversely polarized target, however, ty l s ^ and the 
two effects can be separated. 

The Collins fragmentation function represents a correlation between the quark's trans- 
verse spin and the transverse momentum of the leading hadron formed in the frag- 
mentation process. The Sivers distribution function reflects the correlation between 
the quark's transverse spin and its transverse momentum within the proton. It has 
been shown 114 lL 14211 that both the Collins and the Sivers effects can contribuite to the 
analysing power observed in the Fermilab E704 p j p — >• nx reaction [43]. Very re- 
cently, Am was measured [44] at RHIC at a much higher energy of ^Js = 200 GeV using 
transversely polarized proton beams. The RHIC data could provide new information on 
the Collins and Sivers functions. 



GENERALIZED PARTON DISTRIBUTIONS 



There has been intense theoretical and experimental activities in recent years on the 
subject of Generalized Parton Distribution (GPD). In the Bjorken scaling regime, exclu- 
sive leptoproduction reactions can be factorized into a hard-scattering part calculable in 
QCD, and a non-perturbative part parameterized by the GPDs. The GPD takes into ac- 
count dynamical correlations between partons with different momenta. In addition to the 
dependence on Q 2 and x, the GPD also depends on two more parameters, the skewedness 
£ and the momentum transfer to the baryon, t. Of particular interest is the connection 
between GPD and the nucleon's orbital angular momentum [45]. 

The deeply virtual Compton scattering (DVCS), in which an energetic photon is pro- 
duced in the reaction ep — > epy, is most suitable for studying GPD. Unlike the exclu- 
sive meson productions, DVCS avoids the complication associated with mesons in the 
final state and can be cleanly interpreted in terms of GPDs. An important experimen- 
tal challenge, however, is to separate the relatively rare DVCS events from the abundant 
electromagnetic Bethe-Heitler (BH) background. From the collision of 800 GeV protons 
with 27.5 GeV positrons, both the ZEUS f£(J and the HI O collaborations at DESY 
observed an excess of e + + p — ► e + + y+ p events in a kinematic region where the BH 
cross section is largely suppressed. The excess events were attributed to the DVCS pro- 
cess and the ZEUS collaboration further determined [43] the DVCS cross section over 
the kinematic range 5 < Q 2 < 100 GeV 2 , 40 < W < 140 GeV. Both the W and Q 2 dep- 
ndences of the ZEUS DVCS cross section data are well described by calculations based 
on GPD and on the color-dipole model. 

At lower cm. energies, the HERMES [14811 and the CLAS [49] collaborations observed 
the interference between the DVCS and the BH processes, which manifests itself as 
a pronounced sin0 azimuthal asymmetry correlated with the beam helicity. Another 
observable sensitive to the interference between the DVCS and the BH processes is 
the azimuthal asymmetry between unpolarized e + and e~ beams. In contrast to the 
Beam Spin Asymmetry (BSA) which is sensitive to the imaginary part of the DVCS 
amplitudes, the Beam Charge Asymmetry (BCA) is probing the real part of the DVCS 
amplitudes Analysis of the HERMES e~ data in 98-99 and the e + data in 99-00 
has shown a positive effect for BSA ll5lll . 

QCD factorization was proved to be valid for exclusive meson production with lon- 
gitudinal virtual photons [52]. Such factorization allowed new means to extract the un- 
polarized and polarized GPD. In particular, unpolarized GPDs can be measured with 
exclusive vector meson production, while polarized GPDs can be probed via exclusive 
pseudoscalar meson production. A broad program of DVCS and hard exclusive pro- 
cesses has been proposed ll53ll for the 12 GeV upgrade at JLab. 
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